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With the aim of deriving symmetric hyperbolic free-evolution systems for GR that possess Hamil- 
tonian structure and allow for the popular puncture gauge condition we analyze the hyperbolicity of 
Hamiltonian systems. We develop helpful tools which are applicable to either the first order in time, 
second order in space or the fully second order form of the equations of motion. For toy models we 
find that the Hamiltonian structure can simplify the proof of symmetric hyperbolicity. In GR we 
use a special structure of the principal part to prove symmetric hyperbolicity of a formulation that 
includes gauge conditions which are very similar to the puncture gauge. 
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I. INTRODUCTION 

In the build up to the solution of the binary black hole 
problem in general relativity, a lot of effort was spent 
examining the Einstein equations as a system of partial 
differential equations (PDEs) . Continuous dependence of 
a solution on given data, or mathematical well-posedness, 
of the initial boundary value problem (IBVP) is an es- 
sential property both for meaningful physical models and 
for numerical applications. When GR is decomposed into 
and against spacelike hypersurfaces it becomes a system 
of ten coupled PDEs; six evolution equations and four 
constraint equations. The constraints must be satisfied 
at all times in the development of a spacetime. There are 
few statements about well-posedness for constrained sys- 
tems [l], , but since the constraints are compatible with 
the evolution equations we may consider free-evolution 
schemes, in which for the purposes of analysis the con- 
straints are not assumed; suitable initial and boundary 
data must be given to guarantee their satisfaction. Well- 
posedness of the free IBVP is determined by the charac- 
ter of the principal part of the system and the boundary 
conditions. In particular well-posedness of the IVP is 
guaranteed by strong hyperbolicity. For the IBVP the 
notion of symmetric hyperbolicity may be used to guar- 
antee well-posedness 0, H| ■ 

From the free-evolution PDEs point of view, the gauge 
freedom of GR has two aspects, the freedom to choose 
arbitrary equations of motion for the spacetime coordi- 
nates and the addition of the constraint equations to the 
evolution equations. Both aspects alter the principal part 
and therefore potentially the well-posedness of the IBVP. 
That the well-posedness of the IBVP depends upon the 
choice of coordinates is troublesome, since there may be 
gauge conditions of interest that do not give rise to a 
well-posed IBVP. 

A coordinate choice that has been applied successfully 
in numerical GR is determined by the puncture gauge [a, 
61 . Here we want to derive Hamiltonian formulations that 
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allow for that gauge condition and possess a well-posed 
IBVP as a free-evolution system. Our methods are also 
applicable when other gauge conditions are desired. 

The motivation is twofold. Firstly every Hamiltonian 
system has a conserved quantity. The existence of a con- 
served quantity is also part of the defining property of 
symmetric hyperbolic PDEs. Therefore one may expect 
a connection between the concepts. The relationship be- 
tween symmetric hyperbolicity and Hamiltonian struc- 
ture has not been previously studied by numerical rel- 
ativists. After analyzing several symmetric hyperbolic 
Hamiltonian systems we find that the two energies are in 
general not related. However in some special cases the 
Hamiltonian energy can be used to simplify the analysis 
of symmetric hyperbolicity. 

The second reason to consider Hamiltonian systems are 
their convenient properties. As discussed above Hamil- 
tonian systems possess a conserved energy. Moreover it 
is guaranteed that the time translation map is symplec- 
tic. The literature contains numerical methods that pre- 
serve one of those properties exactly @, [1] , which makes 
Hamiltonian formulations of GR potentially interesting 
for numerical simulation. Previous studies of symplectic 
integrators in numerical GR show that indeed their appli- 
cation can have advantages over standard Runge-Kutta 
schemes 

In section [H] we describe the basic concepts of hyper- 
bolicity for fully second order systems. We then intro- 
duce Hamiltonian systems in section HID In section IIVI 
we present several new tools that are useful in the anal- 
ysis of hyperbolicity. Next those tools are applied to two 
toy models in section [V] We then turn our attention to 
the Einstein equations. Section IVII contains our Hyper- 
bolicity analysis for a large class of Hamiltonian formu- 
lations of GR with live gauges. Finally we conclude in 
section IVII1 

II. HYPERBOLICITY OF SECOND ORDER 
SYSTEMS 

We start by introducing strong and symmetric hyper- 
bolicity for first order in time, second order in space as 
well as fully second order systems. The former case was 
treated comprehensively in 12]. For fully second order 
systems one can find material on strong hyperbolicity 
in [TH, [HI . Here we also discuss symmetric hyperbolicity 
of fully second order systems. 

A. Definitions 

Well-posedness and hyperbolicity: Well-posedness is 
the requirement that an initial (boundary) value problem 
has a solution that is unique and depends continuously 
on the initial (and boundary) data. It can be shown |§, H| 
that the first order evolution system 

dtu" = A^ v d v u v + S" (1) 



has a well-posed initial (boundary) value problem if it is 
strongly (symmetric) hyperbolic. 

The system (fT|) is called strongly hyperbolic if the prin- 
cipal symbol, the contraction with a spatial vector Sj of 
the principal part A pfJ, v s p , has real eigenvalues and a 
complete set of eigenvectors which depend continuously 
on S{. The system is called symmetric hyperbolic if there 
exists a hermitian symmetrizer H such that (suppressing 
some matrix indices) 

HA p = (HA p ^. (2) 

and where H is positive definite. A hermitian H that sat- 
isfies © but is not necessarily positive definite is called 
a candidate symmetrizer. Symmetric hyperbolicity is a 
more strict condition than strong hyperbolicity jl 2l . ap- 
pendix C] . It is equivalent to the existence of a conserved 
positive definite energy 

E = J edx, e = v) Hu (3) 

for the principal part of the system linearized around 
some background. 

These definitions were used to define strong and sym- 
metric hyperbolicity for first order in time, second order 
in space systems [12[ of the form 

d t v = A\diV + Aiv + A 2 w + a, (4a) 

d t w = BfdidjV + B\diV + B-lv + B l 2 ditv + B 2 w + b. 

(4b) 

In what follows we also consider fully second order sys- 
tems of the form 

d*q = A ij d i d j q + B%d t q + S. (5) 

In fact the set of fully second order systems is a special 
case of the first order in time, second order in space sys- 
tems. However in both forms of the equations certain 
calculations are simplified. 

The system @ is called strongly (symmetric) hyper- 
bolic if there exists a reduction to first order that is 
strongly (symmetric) hyperbolic. We likewise call the 
system §5§ strongly (symmetric) hyperbolic if there exists 
a fully first order reduction that is strongly (symmetric) 
hyperbolic. 

With those definitions a system is strongly hyperbolic 
if and only if the principal symbol, P s , has a complete set 
of eigenvectors (with real eigenvalues) that depend con- 
tinuously on a spatial vector Si. Strong hyperbolicity is 
equivalent to the existence of a complete set of character- 
istic variables with real speeds f |12Hl4] | and appendix EJ). 
For first order in time, second order in space systems the 
principal symbol is 

B?s iSj B^ ) ' W 
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and for fully second order systems it becomes 



1 

*/4. ^ s^Sj 13 Si 



(7) 



The existence of a symmetric hyperbolic first order re- 
duction is equivalent to the existence of a conserved pos- 
itive definite energy, E — J edx (again details can be 
found in fl2H"l4l | and appendix O respectively) . For first 
order in time, second order in space systems the energy 
density has the form 



-1 = A H l{ V ) u 3 

<),c N 



djv 



(8) 



and for fully second order systems it is 



e 2 = u 



H v 

diq 

d t q 



(q)uj 
t 



H# H 



22 



d jq 
d t q 



(9) 



In both cases we denote H v a symmetrizer. Furthermore 
if there is a matrix H lJ such that u[H l3 Uj is a conserved 
quantity, but not necessarily positive definite then we call 
H lJ a candidate symmetrizer. 

Since there is only one time coordinate, if a first order 
in time, second order in space system (jlj is obtained as 
the reduction of a fully second order system §5§ , then the 
two systems will have the same level of hyperbolicity. 

In App.|X]we demonstrate the statements about equiv- 
alence of the various flavors of hyperbolicity for fully sec- 
ond order systems. 

We denote the matrices 



n 2 



*"<■' = ( J??,- & ) <*«»> 



the principal part matrices of the first order in time sec- 
ond order in space and the fully second order system 
respectively. 



B. Symmetric hyperbolicity of fully second order 
systems 

Since the structure of the matrices in equa- 

tion ([TO]) is quite simple one can derive another general 
criterion for fully second order systems of the form ^ to 
be symmetric hyperbolic: 

Lemma 1. Symmetric hyperbolicity of the fully second 
order system |5]) is equivalent to the existence of a second 
order symmetrizer and fluxes (Hi, d> 1 , i? 3 ) satisfying 



for every spatial vector Sj, and 
H ij = 



(13) 



HiA ij + B^HiB j) - B(* V) + 0«t ^ - b^Hi 
ft - HiB j Hi 



hermitian positive definite. 

The lemma is proven in appendix [X] To simplify the 
criterion further we notice the following. If T*j is an 



invertible matrix then positivity of (1131) is equivalent to 
T % jH ok T^ being positive definite. We choose 



r£il 



L+ 



(14) 



with Hi = L~ l . This is possible because positivity 
of ([T3"]) implies positivity of H\. The invertible matrix 
can be chosen to be the Cholesky decomposition of 

Hi. 

It follows that positivity of Q13[) is equivalent to posi- 
tivity of 

L- x A l] L + ittfBtL - tfftLtffiL + L^<p ij f L. (15) 

If we redefine <p l & n d 4' lJ appropriately then we get the 
following 

Corollary 1. The fully second order system ^ is sym- 
metric hyperbolic if and only if there exists an invertible 
matrix L and fluxes ((j) 1 , 4> lJ ) satisfying 



(16) 



such that 

SiSjSk^-L^B^ 'L~^)L- x A ik L = 

= l)A }k ^L-\<\> 1 - L^B'^SiSjSk, (17) 

for every spatial vector Si , and the matrix 

L- l A lJ L + <ib {i L- l B j) L - 0(«y> + (18) 

hermitian positive definite. 

The conditions of corollary [TJ simplify significantly in 
the special case that B l = 0. In that case symmetric hy- 
perbolicity is equivalent to the existence of an invertible 
matrix L and fluxes df — dM^ such that 



Hj 1 = L- 1 A ij L + 6 ijii 



(19) 



is hermitian positive definite. We can choose <\> % = in 
corollary [TJ without affecting positivity, because ir- 1 ^' is 
clearly positive semi-definite. 



AJ'it 



(11) 



III. HAMILTONIAN SYSTEMS 



s^s^ - B^Hi)A jh = 

= A^ k (4> i -HiB i )s i s j s k , (12) 



In this article we aim to analyze the hyperbolicity of 
Hamiltonian systems. Here we give a short introduction 
to the basic notions. 
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A. Variational principle 

In many cases the dynamics of a physical system can 
be described by a single functional, the action 



S = / L ^ q '^ >d1, 



(20) 



with the Lagrangian L, where q denotes the time deriva- 
tive of q. 

According to Hamilton's principle classical systems be- 
have such that the action is minimal. This leads to the 
Euler- Lagrange equations 



d dL dL _ (| 
dt dq dq 



(21) 



The ADM equations [15J are of the form fl2TJ, if the ex- 
trinsic curvature, K, is interpreted as an abbreviation 
for the time derivative of the 3-metric, 7, using (|73al) . 
They are composed of the first order in time equations 
((74)l and the second order equations (|73h>[) . The former 
being the constraints and the latter the dynamical ADM 
equations. 

As mentioned in the previous section we are only in- 
terested in free evolution systems of the form ((4]) or 

For general relativity that means we consider systems 
whose solutions can be directly related to solutions of the 
ADM equations if the constraints are satisfied, but that 
also possess constraint violating solutions. 

The term in the Euler-Lagrange equations (|2ip that 
contains the highest order time derivatives is 



d dL 
dt dq 



d 2 L 



d 2 L 



dqdq dqdq 
Hence, for a free evolution system the matrix 

M := 



(22) 



(23) 



is invertible (the equation of motion of every variable 
may be solved its second time derivative). For this type 
of system one can simplify the structure of the Euler- 
Lagrange equations by writing them in their Hamiltonian 
formulation. 



B. Hamiltonian structure 

A Hamiltonian formulation is given by the specification 
of a Hamiltonian % for the system, constructed from so- 
called canonical positions and momenta (q;p). It can be 
shown that the Hamiltonian is a conserved quantity of 
the system, the total energy. 

Here the Hamiltonian can be expressed as the integral 
over space of a Hamiltonian density 



H= d 3 xHD(q,diq,p,dip), 
is, 



(24) 



where the Hamiltonian density is a local function of the 
canonical variables and their spatial derivatives. One ob- 
tains the following canonical equations of motion 



d t q = 
d t p = 



dHp 
dp 
dU D 



0, 



dU D 

d(d iP y 

dH D 



(25a) 
(25b) 



dq d(diq) ' 

A system that can be written in the form (1231) is said to 
have Hamiltonian structure. 

If one treats the Hamiltonian formulation of a model 
as fundamental then also free evolution equations with 
singular matrix M~ l exist. In some situations it may 
even be possible to reduce those systems to a fully second 
order form. We consider only formulations in which one 
can use the standard Legendre transformation to convert 
between (EH) and (l25l). 



C. Canonical equations of motion 

For systems with Hamiltonian structure one may easily 
show that a candidate symmetrizer exists. We consider a 
Hamiltonian formulation of some quasilinear field theory 
in which the functional form of the Hamiltonian density 
is given by 



H D = ]^H(q)u + S(q, diq,p) 



with — (diq',pl), wher 



yij pi t 



(26) 



(27) 



and S(q, diq,p) does not contribute to the principal part 
of the system. Without loss of generality we can assume 
that 

M- Y = M-\ V l ° =V jl \ V lj = V j \ 

The ansatz (1261) might seem very restrictive, because e.g. 
the Hamiltonian density for the dynamical ADM equa- 
tions contains curvature terms, i.e. second order spatial 
derivatives of the 3-metric. However using integration by 
parts, those terms can be transformed to match with the 
ansatz (|26l) . Integration by parts introduces boundary 
terms, but the equations of motion are unaffected. 
The canonical equations of motion are 

d t q = M~ 1 p + F i d l q + s q , (28a) 

d t p = F^dip + V ij djd iq + s p , (28b) 

where s p and s q denote terms that do not contribute to 
the principal part of the system. 

If we identify q and p with v and w respectively then 
we see that (|28p is a first order in time second order in 
space system of the form A short calculation shows 
that the matrix H from (|27[) is a candidate symmetrizer 
for (|28[) . we denote it the canonical candidate. Thus, if H 
is positive definite, like e.g. for the wave equation, then 
the system (|28p is automatically symmetric hyperbolic. 
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IV. TOOLS TO ANALYZE HYPERBOLICITY 

In this section we present several tools and important 
special cases that are helpful in the analysis of hyperbol- 
icity. 



A. Decoupled variables 

The special case that simplifies our analysis the most 
is that of the principal part matrix having a special block 
structure. 

An obvious special case is that of sets (vi,wi) and 
(«2, W2) of fields that decouple in the principal part of the 
system. It is easy to check that the system is strongly 
(symmetric) hyperbolic if and only if the two subsystems 
are. Hence one may analyze two disjoint systems in iso- 
lation. The most important example in the context of 
general relativity is that of matter and spacetime vari- 
ables. 

Another important special case is that where there are 
pairs (vi,w±) and (1)2,1112) that in the principal part are 
coupled through first order spatial derivative terms only 
and where the couplings within each pair occur through 
zcroth and second order spatial derivatives. The princi- 
pal part matrix has the structure 



/ 





AP .3 


^13* 


\ 




AP J 
^211 








^-24 i 




AP 3 










V 





AP j 
^42 


^43 


/ 



(29) 



For many systems that are relevant for numerical relativ- 
ity one can achieve this form, because one can distinguish 
between variables with an odd or even number of indices. 
Since only the metric can contract these indices the two 
groups of variables are usually not mixed. 

It can be checked that if there is a symmetrizer H , 



Tji3 
•"11 

TT3i\ 
n 12 

Tjj\ 

n 13 
tt3\ 



R t3 

Tjij 

tt3\ 
n 23 

rrit 
^24 



H{ 3 H 



'23 
#33 

#34 



#34 
#44 



(30) 



of 



then 



/#ii 








#14 \ 





R ij 
n 21 


#23 








tt'3\ 
n 23 


#33 





V mt 








#44 / 



(31) 



is also a symmetrizer of that system. The positivity of H 
follows, because the positivity of H implies positivity of 
all its principal minors (exchange rows and columns ap- 
propriately). Furthermore it is straightforward to check 
that H defines a conserved quantity, because the princi- 
pal part matrix has the block structure (|29[) . 



Hence a system with a principal part matrix of the 
form (|29l) is symmetric hyperbolic if and only if it admits 
a symmetrizer of the form (131 [) . 

For fully second order formulations the principal part 
matrix is Wi J defined in (jlOp . In that matrix the up- 
per left block always vanishes and the upper right block 
is the identity. Hence, those submatrices automatically 
have the desired block structure (|2T)|) and the argument 
reduces to conditions on the block structure of A tJ and 
B l . If they have the following form 



A ij 



A " 1 % 

^4 on 



B> 



B\ 2 
B\ x 



(32) 



then one can assume that the symmetrizer has the block 
structure ([3"T]) . 

These considerations about the block structure of sym- 
metrizers can be very helpful when one tries to derive a 
symmetrizer by making an ansatz with certain param- 
eters. A block structure like (l3Tj) halves the number of 
parameters in the ansatz, simplifying the calculations sig- 
nificantly. 



B. Searching for symmetrizers in a family of 
matrices — the rank criterion 



In this section we consider whether or not symmetriz- 
ers may be found in a set of matrices Q . We find a nec- 
essary condition for Q to contain a symmetrizer. If Q 
can be chosen big enough then one obtains a necessary 
condition for symmetric hyperbolicity. We refer to the 
condition as the rank criterion. 

We find in applications that it is typically easier to 
check the rank criterion than the alternative, which is 
to take the set of candidate symmetrizers for a certain 
formulation and prove that the set contains no positive 
definite matrix. 

Consider a family of formulations that is parametrized 
by / formulation parameters, c G C l C t 1 . For a for- 
mulation with parameters c we denote the principal part 
matrix A p f '(c) S R kxk (this can be the fully second order 
or the first order in time, second order in space principal 
part matrix (flU)) ). Assume that A p i 3 depends continu- 
ously on c. 

For the criterion to be applicable we also need to know 
that in every neighborhood U c C C l of each point c E C l 
there are formulations that do not possess a symmetrizer 
in the set of matrices, Q. For our examples we find this 
situation, e.g. when the set of strongly hyperbolic for- 
mulations inside C l has lower dimensionality. 

The idea is then to show that by continuity Q can only 
contain symmetrizers for c if the space of candidates for 
c, § c C 5, is bigger than the space of candidates for a 
generic c G U c . 

We assume that Q is linearly parametrized, i.e. there 



exists a linear map G 13 : 



p k X k 



, that assigns her- 



mitian k x fc-matrices to parameters, g £ R n such that 
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Q = G ij (R n ). According to [3 a matrix G ij (g) is a can- 
didate symmetrizer of a formulation c € C', if and only if 
the following matrix is hermitian for every spatial vector 
s 

S l G^(g)AP J k (c)s p S k , (33) 

with 

«.=(?;). (34) 

Since G lJ is a linear map this defines linear equations for 
the parameters g of the form 

B(c)g = 0, (35) 

where B : C l —> ]R mx ™ is a continuous map, because we 
assume that A p i J is continuous. 

In appendix [B] we show the following. If for fixed 
N G N there is in every neighborhood of c a formula- 
tion c with rank(£?(c)) = N such that £7 does not contain 
symmetrizers for c and rank(_B(c)) > A then there is no 
symmetrizer for c in Q . 

If Q is the set of all hermitian k x /c-matrices and the 
rank of B is continuous at c this implies that c is not 
symmetric hyperbolic. Hence, it is not necessary to dis- 
cuss positivity of matrices, one only needs to calculate 
the rank of B. 



C. Symmetric hyperbolicity of special Hamiltonian 

systems 

In section III Bl we derived general criteria for a fully 
second order system to be symmetric hyperbolic. We 
now discuss simplifications that can be achieved when 
the system additionally has a reduction to first order in 
time, second order in space with Hamiltonian structure 
in the principal part. 

Assume that there exist matrices M — M\ F l and 
yii = v&) t suc h that 

d'fq = M-^V^ - F^MF^didjq 

+ [F^r 1 + M^F^Mdidtq + S, (36) 

where S denotes terms that do not contribute to the prin- 
cipal part. 

The underlying Hamiltonian structure guarantees the 
system a canonical candidate symmetrizer 

W 3 =M^ 1^ (37) 

the matrices MB 1 and MA 13 are hermitian for every i,j: 

S? 2 = MB 1 = MF i + F^M, 
tfj = MA 13 = V ij - F {i HlF^ . (38) 



However one needs to be careful. Although the canon- 
ical candidate (|57)l is block diagonal one cannot assume 
that there exists a positive definite symmetrizer with this 
block structure. However if B l vanishes and the system is 
symmetric hyperbolic then it is easy to check that there 
exists a block diagonal symmetrizer (this was discussed 
previously for general fully second order systems). 

There are several special cases where the existence of 
the canonical candidate (I37|) simplifies the construction 
of a positive symmetrizer: 

1. If M and are positive definite, the system is 
automatically symmetric hyperbolic. 

The proof is trivial, because the canonical candidate is 
positive definite. □ 

2. If A lJ , B % and M commute for every then the 
system is symmetric hyperbolic if A lJ has only positive 
eigenvalues. 

M is hermitian and therefore diagonalizable. The ma- 
trices A 10 , B l and M commute, therefore each A lJ and B % 
is diagonal in the basis where M is diagonal. In this ba- 
sis it is then obvious that MT-L is a candidate and if A 13 
has only positive eigenvalues then MT-L is also positive 
definite. □ 

3. The system \3b}) is symmetric hyperbolic if there 
exists a matrix S such that the two matrices 

SM, SMA 13 (39) 

are hermitian positive definite and 

SMB i (40) 

is symmetric. The symmetrizer is then ST-l . 

Positivity of STi. is obvious because of conditions (|3"9")l 
and a straightforward calculation shows that ST-L is also 
a candidate symmetrizer because of P0|) . □ 

V. TOY PROBLEMS 

We now use the criteria of the previous sections to an- 
alyze several toy problems. Our aim is to demonstrate 
the methods without considering equations with a com- 
plicated structure. 

A. Electromagnetism 

The structure of the Maxwell equations is very simi- 
lar to that of the Einstein equations. In particular they 
share a very similar gauge sector. To expand the vacuum 
Maxwell equations in flat-space 

d l B i = 0, d i E i = 0, (41a) 
dtB, =-{dx E)i, d t Ei = (dx B) t . (41b) 

with their largest gauge freedom we introduce the stan- 
dard scalar and vector potentials cj> and A 1 , 

Bi = (dx A) h Ei = -dtAi - d l( j). (42) 
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The gauge freedom is in the time derivative of (j> and the 
divergence of Aj. 

Hamiltonian structure: In the canonical variables 
(A l ;iTi = —Ei) the dynamical part (|41a[) of the system 
has Hamiltonian density 

1 



+ 4>c m , 



(43a) 



and in analogy with GR, the "momentum" constraint is 



= o. 



(44) 



Here 4> is considered a given field. Following the approach 
of [16j with the Maxwell equations, we promote </> to the 
status of an evolved field by introducing a canonical mo- 
mentum 7r. In flat-space the dynamical vacuum Maxwell 
equations with live gauge conditions can be expanded 
and written with Hamiltonian density 



Hd — T~Lm + T~Lg, 
U g = (ci - 1)0C„ 



1 



c 2 tt 2 - c 3 ndiA l 



(45a) 
(45b) 



where Hg denotes the Hamiltonian for the gauge sector 
of the theory. The canonical variables are (A 1 , <f>; iri, ir). 
In principle the c; are arbitrary given scalar functions on 
the spacetime, however from the point of view of hyper- 
bolicity analysis we will treat them as constants. Since 
we are in flat-space we do not denote the densitization 
of the variables. It is straightforward to generalize the 
system to curved space. The equations of motion are 



d t A i = n i - cidV, 
d t <f> = -c 2 tt - c 3 diA l , 
d t TVi = cPdjAi - d i d j A j - c 3 diir, 
d t ir = -ddiir 1 . 



(46a) 
(46b) 
(46c) 
(46d) 



Note that in (|45b[) the appearance of 4>C m results in the 
addition of the momentum constraint to the evolution 
equations and to a rescaling of <fi (we have ()46a[) instead 
of (H2J) with TTi = —Ei). The constraint addition may 
alter the hyperbolicity of the system. Freezing 0, identi- 
fying 7r = r — diA 1 and choosing the remaining constants 
appropriately reduces the system to the KWB formula- 
tion of electromagnetism [171 ] . 

To see that the canonical momentum of <f> should obey 
the constraint 7r = 0, note that we may write 



d t C m = -c 3 did l Tr. 



(47a) 
(47b) 



A natural choice for the free parameters is the Lorentz 
gauge c\ = ca = C3 = 1. In this case equations (|4"6l are 
a first order in time second order in space formulation of 
decoupled wave equations. In what follows we will con- 
sider an arbitrary gauge choice Ci and derive conditions 
that the resulting PDE system be symmetric hyperbolic. 



In the notation of (|?7|) we have u = (diA 1 , di4>, ni, 
and 



' -C2 I ' 



f k 



-cirf" 
-c 3 8 l k 



(48a) 
(48b) 



( (48c) 

Note that the structure of the gauge sector prevents 
the evolution equation of ir from containing terms like 
did l 4> : since the evolution of the constraint subsystem 
must be closed. The index structure of the variables also 
prevents the evolution equation for A % from containing 
terms like did l <fi. These two facts together guarantee the 
empty row and column of V lJ . In fact any system with 
similar structure in the gauge have the same property. 
We will see this for the Einstein equations in the following 
sections. 

The expanded Maxwell equations (|46l) can be written 
in a fully second order form provided that C2 =/= 0- One 
obtains 



d?A l = djd j A i + (cl/c 2 ~ 1) d i d j A j 
- (cic 2 - c 3 ) /c 2 d lj dt(t>, 
= cfedidV + fctca - c 3 )d i d t A i 
and we read off the principal part 

A ijl k= ( V ij 5 l k + (i-l)v l{i 6 j h 



c\cir\ %] 





(cic 2 - c 3 )S z k 







(49a) 
(49b) 



(50a) 
(50b) 



Strong hyperbolicity: Following [1 21 ] we investigate 
strong hyperbolicity by performing a 2 + 1 decomposition 
in space. The principal symbol then decomposes into a 
scalar and a vector block. The characteristic speeds in 
the vector block are ±1 and this submatrix is always diag- 
onalizable. In the scalar block of the expanded Maxwell 
equations the characteristic speeds are ±^cic 3 . Thus, 
strong hyperbolicity requires that c\ and c 3 share a sign, 
or one of them vanishes. Further calculations reveal that 
the system is not strongly hyperbolic in the latter case. 

Strong hyperbolicity additionally demands a complete 
set of characteristic variables. They exist only if cic 2 = 
c 3 7^ 0, i.e. if in the fully second order system (|49[) the 
vector and scalar potentials are decoupled. The charac- 
teristic variables are given by 



U s±Vl = d t 4> ± ciyfc^dscf), 
U s ±v 2 = d t A s ± c ly /c2~d a A a 
U A ±i=d s A A ±d t A A . 



(51a) 
(51b) 
(51c) 



The non-trivial characteristic speeds are V\ — i> 2 = 
ici J02 so strong hyperbolicity requires c 2 > as well. 
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Symmetric hyperbolicity: Before we start the con- 
struction of positive symmetrizers we check whether the 
necessary conditions for the existence of positive candi- 
dates, i.e. strong hyperbolicity and the rank criterion, 
are satisfied. Thus, we immediately demand strong hy- 
perbolicity so B l vanishes, C2 > 0, C3 = c\Ci ^ 0. Appli- 
cation of the rank criterion (see section UV Bp results in a 
similar restriction. We find that the rank of the relevant 
matrix B{c) changes if and only if C3 = c\Ci. 

The remaining matrix A l ° is automatically symmetric. 
So, taking Hi = 1, (j> % — 0, y = in lemma [TJ the sys- 
tem is obviously symmetric hyperbolic if A lJ has positive 
eigenvalues. Positivity of the eigenvalues with C1C2 = C3 
is equivalent to the condition \ < c\c 2 < 3. In geomet- 
ric units the characteristic speeds v of the system satisfy 
l/y/2 < \v\ < s/3. 

Instead of the canonical candidate we may also con- 
sider the case where 



h m h = ( hn 







= W, 



(52) 



with some constant b\ . The eigenvalues of the resulting 
candidate symmetrizer are 



Ai 



1 + c\c 2 ~ h 



A, 



3 - c\c 2 + h 



2 2 



A 3 = -1 + 2c x c 2 + bi, A 4 =c 1 c 2 , 
A 5 = 1. 



(53) 



Hence, if c 2 > we can always choose b\ such that 
this candidate is positive definite. If c\c 2 < 4/3 then 
b\ = 1 is a possible value and for c\c 2 > 4/3 one can 
take bi = c\c 2 . Thus, every strongly hyperbolic formula- 
tion of Electromagnetism with Hamiltonian structure is 
symmetric hyperbolic. 



B. The pure gauge system 

In this section we consider scalar fields (T, X 1 ) with 
second order equations of motion. We determine the set 
of equations of motion with Hamiltonian structure which 
are strongly and symmetric hyperbolic. Finally we con- 
sider which equations of motion obtained for the lapse 
and shift when the fields are taken as coordinates. 

Hamiltonian structure: We start by introducing 
canonical momenta (0,Hi) for the scalar fields (T, X 1 ) 
and write the most general ansatz for the Hamiltonian 
without introducing new geometric objects such that the 
equations of motion are linear and naturally written with 
"9o = h(dt — P l di) v derivatives (where a and j3 l are lapse 
and shift of the background manifold slicing respectively) 



2^/j 2 2^f 

+ ^y/lD i XW i X j + c J^^ Dl XW 3 X l 



(54) 



CeaOD.X 1 + c 7 aH l D l T + Qp i D i T + li^D.X 1 



Terms like QX l DiT would make the equations of mo- 
tion nonlinear. We could also adjust the system by ad- 
ditional source terms. In the notation of ((27)) we have 
u = (djX k ,d t T,ni,Qy and 



M 



-Ilk 



a 



, Ik 



C37""' 
2x/7 V ci 



r k — — 



C67 



c 7 8 l k 
a o 



/3 l ( 5' k 
1 



V 



2 V c 27 i ' 



kl 



(55a) 
(55b) 
(55c) 



The pure gauge system does not carry any constraints, 
so unlike in the Maxwell equations the V lJ u matrix may 
contain terms on both diagonal components. Otherwise 
the block structure guarantees that the the matrices look 
very similar to those of the Maxwell theory ([4"5]) . In fact 
electromagnetism is a special case of the pure gauge sys- 
tem. 

Fully second order system: The number of parame- 
ters now present in the system complicates the analysis 
of hyperbolicity. To simplify the calculations we imme- 
diately switch to the fully second order version of the 
system, in which the equations of motion have principal 
matrices 



B l = aB l n 



2/3 1 



(56a) 
(56b) 



where the principal part in normal slicing (a = 1, j3 i = 0) 
becomes 



C17 



C4<^ 
CB**J 



Ao^k = 

B l ji = 
and we write 

Cl = C3C4, 

C3 = — (c 2 c 3 - c 7 ), 

C3 
1 , 

C5 = — (ClC 7 + C 3 C 6 ), 

C3 



C27' 









C37* 



(57a) 
(57b) 



C2 = — (cic 5 - c%), (58) 

^1 



C4 



1 

Cl 



(cic 7 + c 3 c 6 ), (59) 
(60) 



to simplify the matrices. In what follows we assume nor- 
mal slicing, i.e. a — 1, (3 l = 0, Ao = A and Bo = B. 

Since we require M~ l invertible, ci and C3 must be 
nonzero, which means that if one of C4 and 65 vanish, 
they both must. 

Given the linear fully second order system (|57[) one 
can always derive a Hamiltonian formulation, provided 
c 4 ^ ^ c 5 , using 



ci 



C4 



C3C5 
c 4 

£1 

C3 



C2 



C3C4 + c 5 c 7 



C3C5 



c 6 = c 5 1 



£7 

C4 
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and 



C5 



C3C4 



C 2 C 4 



c\c 5 



2c 4 c 5 c 7 + C5C7) , 



together with (|55[) where C3 and C7 can be chosen freely. 

Strong hyperbolicity: We assume a generic choice of 
the parameters ci,...,5s, i.e. no degeneracies in the 
characteristic variables occur. Then we perform a 2 + 1 
decomposition in space and introduce the auxiliary quan- 
tities 



Ax = VC1+C2, A T = y/c3, 



(61) 



the speeds of the decoupled wave equations (in the scalar 
sector) in normal slicing when B vanishes. In the vector 
sector the fully second order characteristic variables are 



Ua±x = d t X A ± \ v d s X A . 



(62) 



with speeds Ay = \fc[, so c\ > is required. With 
non-vanishing B the characteristic speeds in the scalar 
sector ±(vt,vx) are modified to 

2v\ = A T + \ 2 X + C4C5 

+ \J (A T - A x ) 2 + c\c\ + 2c4C 5 (\ 2 x + A T ), (63a) 



2v x = A^, + A 



x 



C4C5 



- \A A ^ " X x) 2 + 4% + ^c 5 (X x + A T ). (63b) 

Note that if C4 or C5 vanish we recover the decoupled 
speeds. The characteristic variables in the scalar sector 
are 



U s ±v x = d t X s ± v x d s X s 



± 



C5VX 



(d t T±v x d s T), 



U s±VT =d t T±v T d s T 
C4VT 



± 



*x 



(d t X s ± v T d s X s 



(64a) 
(64b) 



when v x =/= At and vt 7^ Ax- Further calculation using 
the fact that v\v\ — X\\\ reveals that v x — At if 
and only if vt = Ax, which is only possible if C4 and 
C5 vanish. In this case the variables decouple and the 
only condition for strong hyperbolicity is that the speeds 
are real. In summary, the system is strongly hyperbolic 
provided that 



\ 2 V = Cl > 0, 
= c 3 + 0, 
v x >0. 



\ X =C 1 +C2^ 0, 

v\ > 0, 



(65) 



To obtain the characteristic variables for the fully second 
order system in a general slicing from the expressions in 
normal slicing one can e.g. simply replace the dt deriva- 
tive in those expressions by a "do" derivative: 



d t X ->■ - (d t X - (5 s d s X) . 
a 



(66) 



E.g. in the vector sector one obtains 



U A 



- (d t X A 
a 



{(3 s T a\)d s X A ) . 



(67) 



Given a characteristic speed v in normal slicing it be- 
comes av + (3 s for general slicings. 

Symmetric hyperbolicity: As in the case of electro- 
magnetism we first demand strong hyperbolicity. The 
second necessary condition, the rank criterion, cannot be 
applied here, because in the parameter space we do not 
find a dense set of formulations which are not symmetric 
hyperbolic. 

We ask whether the modified canonical candidate 



fei 



M M 



with the matrix </> y ki 
immediately get that 

c 3 > 0, 



(68) 



defined in (|52|) is positive. We 



c 4 /c 5 > 0, 



c 3 > (69) 



must be satisfied, because the candidate is block diagonal 
with one diagonal block which has those elements. 
The eigenvalues of the remaining upper left block are 

2ci 2ci + c 2 , ci + 2c 2 , „, ^ 



01 



2bx. 



2c 3 2c 3 c 3 

Hence, we can choose b± such that the modified candidate 
becomes positive provided that 



ci > 0, 



ci + c 2 > 0. 



(71) 



Strong hyperbolicity of the system only implies that ci > 
0, A T A X > and (A T - A^) 2 + cl^ + 2c 4 c 5 (A| + A T ) > 0. 
Hence, as opposed to Electromagnetism, we get strongly 
hyperbolic formulations where the modified canonical 
candidate is not positive definite. 

However, here the expressions are still simple enough 
to consider all matrices that can be constructed from 
the metric. We derive general candidate symmetrizers in 
that class and analyze their positivity. We find that for 
a generic choice of the parameters every strongly hyper- 
bolic formulation is also symmetric hyperbolic. 

Discussion: It is not clear by examining matrices (|57p 
exactly what restriction Hamiltonian structure puts upon 
the pure gauge system. The only condition is that when 
either of C4 or C5 vanish then they both do. One may 
analyze the level of hyperbolicity of the system even if 
Hamiltonian structure is abandoned. The two cases C4 = 
and C5 = are similar, so we describe only the first. We 
find that whenever Ax 7^ At the system is both strongly 
and symmetric hyperbolic whenever At, Ax > 0. 

With this example we demonstrate that there are 
strongly hyperbolic Hamiltonian formulations where ap- 
parently the existence of the canonical candidate does not 
simplify the symmetric hyperbolicity analysis. Moreover, 
in the fully second order form of the equations of motion 
Hamiltonian structure does not play a role. Thus, if gen- 
eral relativity is considered then we do not expect to 
find symmetrizers through a simple modification of the 
canonical candidate. 
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Evolution of the lapse and shift: To derive equations 
of motion for the lapse and shift from our pure gauge 
system we may simply use the relationship between the 
lapse and shift a, f3 l and the local coordinates in time 
and space (t, x z ) 

a = —±r-, P l = -an a d a x\ (72) 
n a d a t 

with unit normal to the hypersurface n a , and take T = t, 
X % = x 1 . The lapse and shift evolution equations do not 
inherit the same principal part as the pure gauge system 
from which they were derived. Furthermore, it is not pos- 
sible to give the fully second order equations of motion 
for the lapse and shift until the gauge is coupled to par- 
ticular equations of motion for the metric and extrinsic 
curvature. 



where the constraints must be rewritten in terms of the 
canonical ADM variables. 

The constraint evolution system is closed, so in the ab- 
sence of a spatial boundary they need only be satisfied 
in one time-slice to guarantee their satisfaction later in 
time. This justifies analyzing the PDE properties of the 
system without assuming the constraints. When a spa- 
tial boundary is present suitable constraint preserving 
boundary conditions are also required. 

In the ADM formulation the lapse a and shift f3 l are 
thought of as given fields. The formulation is weakly, but 
not strongly hyperbolic (the principal symbol has real 
eigenvalues) . The initial value problem is not well-posed. 



B. Gauge conditions 



VI. FORMULATIONS IN NUMERICAL GR 

In this section we investigate formulations of GR with 
Hamiltonian structure. The aim is to derive a symmetric 
hyperbolic Hamiltonian formulation with properties that 
are believed to be important for stable numerical simula- 
tions and to relate this formulation to currently existing 
ones. 

A. The ADM system 

The starting point in the derivation of a new formula- 
tion is the ADM system [l5[ . It can be derived through 
a 3 + 1 decomposition of the vacuum Einstein equa- 
tions 18]. If one writes the result in geometric variables 
(jij , Kij ) , where 7-y is the 3-metric and Kij is the ex- 
trinsic curvature induced in a spacelike slice, then one 
obtains equations of motion 

dtli] = -2aKij + Cpjij, (73a) 
d t Kij = -D t D 3 a + a[Rij - 2K lk K) + K t] K] 

+ C p K tJ . (73b) 

Equation (|73al) defines Kij and eqn. (|73b[) is the projec- 
tion of the Einstein equations out of the slice. By con- 
tracting into the slice one also obtains the Hamiltonian 
C and momentum C, constraints, which are given by 

C = R + K 2 - K tJ K l] = 0, (74a) 
d = Dj (K ij - j ij K) = 0. (74b) 

The ADM equations can also be written in terms of 
canonical ADM variables, (7^ ; ir^), through the relation 

ir lj = y^f(KY j - K ij ). (75) 

In these variables the dynamical ADM equations ((73)) 
have a Hamiltonian structure with Hamiltonian 

Hadm= f(-aC + 2(3 l C t )^d 3 x, (76) 



In the introduction we noted that from the free- 
evolution PDEs point of view, the gauge freedom of GR 
has two aspects, the first of which is the ability to choose 
arbitrary equations of motion for the lapse and shift. 
Choices in which either the lapse or shift evolve dynam- 
ically arc known as live gauge conditions. In this section 
we describe two popular live gauge conditions. 

Harmonic gauge: In terms of the lapse and shift the 
harmonic gauge condition DX° = becomes 

d t a = -a 2 K + /3 p a p , (77a) 
dtft = Pdrf 1 - ad 1 a + a 2 j jk r jk . (77b) 

The use of this gauge was the key component in the first 
proofs that GR admits a well-posed initial value prob- 
lem [Tj|, since when it is appropriately coupled to the 
EE's they become a particularly simple symmetric hyper- 
bolic system of wave equations. When the gauge is mod- 
ified with arbitrary source terms hyperbolicity of that 
formulation is not affected [13] • The modified gauge is 
known as generalized harmonic gauge or GHG. 

Puncture gauge: We call the combination of the pop- 
ular Bona-Masso lapse and gamma driver shift conditions 

d t a = -Li L a 2 K + f3 p a p , (78a) 

dtfr =n s r i + P j d j f3 i -r l i3. (78b) 

with 

fi = 7 i/3yfe r ; fc + l~fi/^T k kj (79) 

the puncture gauge. In fact there are various choices of 
parameters inside the conditions which we do not con- 
sider here; for a review of the choices made by vari- 
ous groups see table 1 of Our gamma driver condi- 
tion (|78bp can be obtained by integrating the most com- 
monly implemented condition once in time Q. Station- 
ary data for the lapse condition was discussed in pll - [23| . 
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C. Hamiltonian formulations 

Now we are interested in whether or not strongly or 
symmetric hyperbolic Hamiltonian formulations with a 
large class of live gauge conditions can be derived. Our 
Hamiltonian formulations are described in terms of the 
canonical variables (7^, a, /S 1 ;^ ,<r, pi). For solutions of 
the Einstein equations a and pi are constrained to vanish. 

To get a Hamiltonian formulation in other variables, 
e.g. densitized ones, one can apply a symplectic trans- 
formation to obtain the correct Hamiltonian. This does 
not alter hyperbolicity. In general, if the desired position 
variables are f(q) then the corresponding canonical mo- 
menta are f'(q)~^p, where /' is the Jacobi matrix of / 
(/ must not contain derivatives of q). 

To derive formulations we follow the approach of 
Brown [l||, i.e. we add appropriate gauge terms to the 
Super- Hamiltonian "Hadm: 

U = Hadm + n GHG + J d 3 x (Aa + W Pi ) , (80) 

where Hadm + T~Lqhg is the Hamiltonian for Brown's 
generalized harmonic formulation jlfjj. The terms that 
appear in the principal part are 

n GHG = FaDia + ftpidifi + a 2 Pl T) k ^ k - ap l D,a 

+ ( - 4a 2 7y W + a 3 a 2 - 4a 3 Pl p^) . (81) 
°V7 

The terms in A and fl l that affect the principal part of the 
equations of motion are linear in the canonical momenta 
and in the first spatial derivatives of the positions. 

I 



lvl klmn 



ij kl ran 



where 

yijklmn _ _ a ^ij^kl^mn + y(fcy)(myi)j 
_|_ i(m^n)(k l)j _ ^ij^k(m^n)l 

Various entries of V % i kl mn vanish for exactly the same 
reason as those in the Maxwell system (|4"5)l . 



We restrict the possible terms further by considering 
those gauges for which the shift appears in the principal 
part only through its spatial derivatives and the obvious 
advection terms. 

This is of course a restriction of the admissible gauges. 
But the live gauge conditions that are used in numerical 
GR are usually of this form. Moreover in those gauges the 
hyperbolicity analysis is simpler, because, as we will see, 
one may linearize around Minkowski space without loss 
of generality. Furthermore the discussion of Sect. lIV A"l is 
applicable. 

One also finds strongly hyperbolic formulations with- 
out this "<9o" restriction, but they are not considered 
here. We get the following expressions 

A = -Cia 2 7 - 1/2 7^ + CWV 1/2 <7+ 

+ C 7 (aD^ - ^aj jk ^da jk ) (82) 

and 

n l = C 2 a 2 T) kl 3k + C 3 a 2 T k k]1 ^+ (83) 
- C^fWja - Ce-f- 1/2 a 3 p 3 -f lj . 

We denote Co , . . . , C7 formulation parameters and ob- 
tain Brown's Hamiltonian [l6j when C7 = and we re- 
place 

Ci^d-1/2, CWCa-l, C 4 ^C 4 -l/8, 
C 5 -> C 5 - 1, C 6 -> C 6 - 1/2. (84) 

For those formulations we get in the notation of (|2"T)l 
that u = {d l ^ km ,d l a, di/3 k , ir km , a, p k )^ and 



(85a) 



(85b) 



(85c) 

I 

To see that for the hyperbolicity analysis one may 
linearize around Minkowski space observe that when 
A p (det(jij),a, f3 l ) is the first order in time, second or- 
der in space principal part for an arbitrary background 
solution and 

T = diag(a- 1 /2 7 -i/4 5a i/2 7 -i/4 ;a i/2 7 -i/4 j 

a- 1 /2 7 l/4 5 a -3/2 7 l/4, a -3/2 7 l/4) (gg) 



(27fc(m7n)/ ~ Jkl7mn -(1/2 + Ci)a~{ kl 
-(l/2 + d)a 7mn (l/4 + 2C 4 )a 2 


P&pltf 28{ llk)n 

pi CjaSi, 

a 2 ((l + C 2 )7 j(m 7 n)fc -5(C , 2-C , 3 + l)7 J ' £ 7 mn ) ~(l + C 5 )a^ k pS k n 

7 1 / 2 I ji( m ^n)j _ ryij^rnn g q 
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then 

r- 1 A p (det(7y),a,/0 i )T = aA p (l, 1,0) + (3 P 1. (87) 

Diagonalizability and eigenvalues are not altered by a 
similarity transformation, i.e. to analyze strong hyper- 
bolicity one may consider the matrix s p A p (l, 1, 0) for an 
arbitrary spatial vector s. 

Concerning symmetric hyperbolicity a short calcula- 
tion shows that when H is a symmetrizer of A p (l, 1, 0) 
then T^HT is a symmetrizer of A p (det(7y) , a, f3 l ) for an 
arbitrary background metric. In what follows we assume 
a = 1 = 7 and /3* = 0. 

Strong hyperbolicity: Again we perform a 2 + 1 de- 
composition in space (l^ |. and the principal symbol de- 
composes into three submatrices, a scalar, vector and 
trace-free tensor block. The characteristic speeds in the 
trace-free tensor and vector sector are 



\ 2 



1, Ay = 1 + C 2 

respectively. In the scalar block one finds 

2v% = X\+X 2 _- C 5 C 7 



(89a) 



+ yj (X\ - X 2 _Y + C 2 C 2 - 2C 5 C 7 (X 2 + + XI), 



2vt 



where 



— A + + A_ — C5C7 



(89b) 



- ^(A* -X 2 _) 2 + C 2 C 2 - 2C 5 C 7 (X 2 + + X 2 -), 



2X 2 ± = (2Ci + C 2 + C 3 - C 7 + 2) 

±(-2C 1 + C 2 + C 3 + C 7 ). (90) 

Note that the speeds in the scalar sector are very similar 
to those of the pure gauge system, and simplify signifi- 
cantly in the special case C5 = or C 7 = (or both). 
They simplify further when A + = A_ . 

Using the techniques described in App. [Owe can clas- 
sify the formulations with respect to their strong hyper- 
bolicity. We find three families ^=1,2,3, of strongly hy- 
perbolic fully second order formulations. 

From the vector block of the principal symbol it is 
obvious that the condition 

C 6 = (Xy - l)/2 (91) 

must be satisfied in every case. Moreover it is clear that 
the characteristic speeds must be real, i.e. 



Xl>0, 



v +/- 



> 0. 



(92) 



The three families correspond to different solutions of the 
requirement that the scalar block is diagonalizable. 

Family T\ : The first family has four free parameters, 
(A+, A_, C5, C 7 ), and one obtains strongly hyperbolic for- 
mulations if A_ 7^ 0, v 2 7^ v 2 , 

0^(1 + 2C 5 )X 2 _ - X 2 + , (93a) 

\2 \2 C^X 2 ^ + 3C 7 X 2 ^_ 

2A y = 3A + - 1 + n ■ or ,_ u2 _ to C 5 , 



8C 4 = A - 1 



(1 + 2C 5 )A^-A2 
(1 + 2C 5 )X 2 _ - A 2 . 



C 7 . 



(93b) 
(93c) 



where we have replaced C±, C2 and C3 using the defini- 
tions jnOD and (JHHJ) of A± and Ay. 

Family Ti : If the inequality (|93a[) is not satisfied, but 
A_ 7^ 0, v\ ^ v 2 _ still hold, then the solution be- 
comes singular. One obtains another family with three 
free parameters, (A_, Ay, C 7 ): 

12C 4 = 3A^ - 3 - 9C 7 - (1 + 6C 7 )(Ay - 1), (94a) 
3C7 . o A 2 



C 5 = - 



1 + 6C7 : 



Ai = 



It must satisfy 



C 7 ^0, 



1 + 6C 7 ' 



C 7 + -1/6. 



(94b) 



(94c) 



The first of these inequalities guarantees that there are 
no more repeated speeds and the second prevents the so- 
lution from becoming singular, both of which are special 
cases that contain no further strongly hyperbolic formu- 
lations. The characteristic speeds in the scalar sector 
become 



2(1 + 6C 7 )4 = 3C 7 (2A 2 _ + C 7 ) + 2A 2 



±C 7 J9(2X 2 _ +C 7 ) 2 + 12X 2 _ . (95) 



Family F 3 : If the scalar block has only one eigenvalue 
(ii 2 . = v 2 ) one obtains a one parameter family of formu- 
lations with A+ > 0: 



A- = A+, C 5 = 0, C 7 = 0, 

\- — •> \- 1 r< _ 1 /vi \- 1 , (96) 



2X 2 v = 3X 2 + -l, C 4 = l/8(X 2 + - 1). 



This is the limit of T\ ([93]) for A + — > A_ along curves with 
C5 = = C 7 . But if one takes in ([53"]) the limit v + —> v_ 
then the result is in general not strongly hyperbolic. To 
obtain ([TO]) one may also take the limit C 7 — > 0, Ay — >■ 
(3A 2 . - l)/2 of F 2 ([Ml) (in any order). Again the case 
C 7 = is not strongly hyperbolic for general choices of 

Ay. 

Fully second order system: The fully second order 
equations of motion can be constructed for each of three 
strongly hyperbolic systems JFi. We present only T 3 . In 
terms of variables u^i = (7fcz, a, Pk) the two parameter 
strongly hyperbolic system J-3 has the fully second order 
principal part 



8 2 u kl = A l] kl mn d t d u mn + B\r n d t d t u mn . 

The principal matrices A^ ki mn and B % ki mn are 



(97) 
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A\\ k i mn <7 B 7 tf 7w \ / 

A ij k i mn = ( C 67 « , S 4 fci m ™ = | , (98a) 

C 7 Y j S^ + C s j m ^5§ / V B 2 j ik j mn S37* 



where 

A\ j lk i mn = Y ] S^ + Ci7 i(m 7" )j 7fci 

+ <V J 7fca m " + <?3<^ )7 mn 

The parameters in A % l 1 hi mn are given by 

2C 1 = 1-A 2 + , 2C 2 = A 2 + -l + 2 ^ + ^ 2 + A 2 M 

(100a) 

2C 3 = 1-A 2 h , 4C* 4 = 3(A 2 ;-1). (100b) 
The remaining parameters in A lJ ki mn are 
Ai -A 2 . 



C 5 = 2- + 



1 + 3A 

4CV = 3A^ + 1, 4f 



._> 6*6 — A_, 

A 2 +1 



(101a) 



(101b) 



The parameters of B l M ran are 



£1 = 2 

B 3 = 3 



A^-Aj 
1 + 3A| ' 

A^-A 2 , 
1 + 3A 2 ' 



A^ - A 2 
1 + 3A 



2£ 2 



(102a) 
(102b) 



Characteristic variables: The characteristic variables 
may be constructed for each family T\. For brevity we 
present them only for the subset C5 = C7 = of the 
first family T\ . The characteristic variables in the scalar 
sector are 

U Sa ±\- = (ft ± \-d s )a =f 2A_(ft ± X-d s )j ss 

T i\-{d t ± \-d s )~/ qq , (103a) 

U S p±\- = (ft ± \~d s )a =F 2A_(ft ± A_ft) 7ss 

T 4A_(ft ±A_ft) 7gg , (103b) 

2U S0±X+ = 2(ft ± A+ft) 7ss T A+(ft ± A+ft)/3 s , 

(103c) 

^s±a + = (ft ± A+ft)( 7ss - 7«)» ( 103d ) 

with speeds A+ and A_ as given above (|89|) . The vector 
sector has 

Ua±j = (ft ± Ayft) 7As , (104a) 

C/a±/3 = (ft ± Xvd s )/3 A . (104b) 



with Ay = \/l + C*2- Finally the tensor sector has char- 
acteristic variables 



C/JI ±1 = (ft±ft)7^ 



(105) 



where 7JJ3 denotes the transverse trace free part of the 
metric against Sj. 

To find the characteristic variables of the system in an 
arbitrary background, one must transform the eigenvec- 
tors of the principal symbol through the similarity trans- 
formation defined by T from (|86[) and replace the time 
derivatives ft — ► ft. 

Symmetric hyperbolicity: As for the Maxwell equa- 
tions we impose strong hyperbolicity and check where 
the rank criterion is satisfied. Interestingly we find again 
that the rank of the relevant matrix B(c) from section 
IIVBI is discontinuous (and drops) at every strongly hy- 
perbolic formulation. Hence, this does not lead to imme- 
diate restrictions in the parameter space. 

A complete analysis of symmetric hyperbolicity seems 
to be too hard in the general case, because expressions 
become very complicated. We want to discuss special 
cases where this problem can be avoided and conjecture 
about the general case. 

In Add. [Dlwe construct positive definite symmetrizers 
for formulations whose parameters satisfy 



A 



v 
C 6 



(3A 2 



3(A 2 



l)/4, 
l)/8, 



C 7 



(A 2 - 
0. 



l)/8, C 5 = 0, 



Ai > 0, 



Ai > 0. 



(106) 



Hence, we found a two parameter family of symmetric 
hyperbolic formulations. Using the same techniques we 
are also able to prove symmetric hyperbolicity of other 
two parameter families. With (I106P they share that they 
belong to the four parameter family T\ (f9"3"| and that 
C 5 = 0. 

For the three parameter subfamily of T\ that satisfies 
also C5 = we can construct candidates such that 
the contractions SiH l3 Sj are positive definite for every 
spatial vector s. We think that the full matrices if y can 
be chosen positive definite in this case, but are not able 
to prove it. 

We also deal with the question whether also in GR ev- 
ery strongly hyperbolic formulation is symmetric hyper- 
bolic, as it was the case for the toy problems in section fVl 
We find that the (strongly hyperbolic) two parameter 
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\ 2 + = (l + C 5 )A 2 



subfamily family of J-\ with 
A 2 , = i(l + C 5 )(l + 3A 2 _), 

C 6 = ^(C 5 (3X 2 _ + 1) + 3(A 2 _ - 1)), C 7 = 0, 

A^ > 0, C* 5 > -1 (107) 

and a generic choice of A"! and C5 is not symmetric hy- 
perbolic. 

To prove that we first search for candidates in the set 
of matrices that can be written in terms of the metric. 
We find a six dimensional space, Q = Gs ®Ga, which de- 
composes into a four dimensional Gs, the symmetric part 
5™', and a two dimensional Ga, the antisymmetric (f>^ 
(|TT|) . We can show that there is no candidate with posi- 
tive definite principal minor SiH^Sj — SiH^Sj (where 
s is any spatial vector). That means the four degrees of 
freedom in Gs are n ot sufficient to construct a positive 
symmetrizer for the principal symbol. Thus, since all 
principal minors of a positive definite matrix are positive 
as well, there is no symmetrizer which can be written in 
terms of the metric. 

We then ask for the candidate symmetrizers in the 
most general set of matrices, i.e. the set which is 
only restricted by the appropriate block structure (1311) . 
There we find a 76 dimensional space of candidates, 
G = Gs © Ga, but again Gs is four dimensional. The 
most general set of matrices of course contains all ma- 
trices that can be written in terms of the metric, i.e. 
Gs = Gs- Hence, G does not contain a matrix with posi- 
tive definite principal minor SiH^ Sj, and therefore there 
is no positive definite symmetrizer for the formulations 
P71> . 

This statement holds for a generic choice of the param- 
eters in p07[) . To deal with special cases we apply the 
rank criterion. It tells us that within (|107[) one cannot 
find symmetrizers, except when Gs is at least five dimen- 
sional. This is the case for C5 = only, which results in 
the (symmetric hyperbolic) one parameter family J-3. 

Thus, the obvious conjecture is that formulations 
where the dimension of Gs is at least five are symmetric 
hyperbolic, and that formulations where Gs is at most 
four dimensional are not. 

To prove this one would need to analyze positivity of 
candidates for three and four parameter families. A prob- 
lem that we could not solve so far. 

We find that the dimension of Gs is at least five e.g. 
when one imposes C5 = in the family T\ (|§3"]l . also for 
Cj 0. But, as discussed above, for this three parameter 
family we can only show that certain sub blocks of a 
candidate are positive definite. 

Discussion: We have also constructed the character- 
istic variables for the four T\ and three parameter T2 
families of strongly hyperbolic Hamiltonian formulations 
of GR. It is interesting that in every fully second order 



strongly hyperbolic Hamiltonian formulation (Maxwell, 
pure gauge and GR) we find that the fully second order 
characteristic variables always take the simple form 



(108) 



with speeds ±i> given constants ay and primitive vari- 
ables qi. There are strongly hyperbolic fully second or- 
der systems without this property, for example the Z4 
formulation coupled to the puncture gauge [2J]. At the 
moment it is not clear what causes the special form (|108l) . 

We originally aimed to find strongly and symmetric 
hyperbolic Hamiltonian formulations with popular gauge 
conditions. Having analyzed a large class of gauges, we 
are now in a position to return to that question. The 
generalized harmonic gauge is straightforward and as dis- 
cussed previously is recovered with Cj = 0, i = 1, . . . , 7. 
More interesting is the puncture gauge (1751 . 

The formulation parameters in the strongly and sym- 
metric hyperbolic formulations can be used to adjust the 
characteristic speeds. If we regard them as scalar func- 
tions of the position variables (7^, a,/? 1 ) then the equa- 
tions of motion will contain derivatives of the parame- 
ters, but in the principal part only the functions appear. 
Hence, hyperbolicity of the formulations is not altered. 

Here we want to present the evolution equations for 
lapse and shift for a formulation that is very close to the 
puncture gauge. It belongs to the symmetric hyperbolic 
family f|106[) . We choose 



4/! S 7 



1/3 



3a 2 



Thus, symmetric hyperbolicity requires 



Ml > 0, 



47 1/3 A*s > a 2 



(109) 



(110) 



which reduces to < a < a/37 1 ^ 6 for the popular choice 
fiL = 2/a, /is = 3/4. To make the comparison with 
other systems simpler we present the equations in terms 
of the Z4 variables. Instead of the canonical momenta 
(ir v a, pi ) the Z4 formulation [25| (which is not Hamil- 
tonian) takes the extrinsic curvature and two fields called 
(0, Zi) as evolved variables. The mapping to the canon- 
ical variables is defined through d75t and 



e = i— 

2^7' 



2^7' 



(111) 



In these variables, restricting to Brown's Hamiltonian re- 
sults in a system that is as close to the Z4 formulation 
as possible. With the choice (|109l) we find 

d t a = (3 l d ia - n L a 2 K + ^ L a 2 @, 

d t p = pd;? + M57 1/3 r;- fe 7 Jfc + kvsi 1/3 - « 2 )rV 3 



+ 2psl 1/3 Z l - aD l a. 



(112) 
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Hence, near the puncture (for a — > 0) we obtain the 
puncture gauge condition ([75]l when the constraints = 
and Zi = are satisfied. 

Brown finds in [l6j that a strongly hyperbolic Hamilto- 
nian formulation with the puncture gauge does not exist. 
Here we have shown that with a slight modification of 
the T-driver shift condition even symmetric hyperbolic- 
ity can be achieved. The new terms are small near the 
puncture. 

The final measure of how useful the formulations with 
these gauges are will be determined by numerical exper- 
iments. Some of the ingredients for successful numer- 
ical evolutions are now understood. In this work we 
have taken care of questions related to the continuum 
formulation. For long-term stable numerical evolutions 
of puncture data a convenient choice of variables is an 
important feature as well [26j]. Furthermore, even if con- 
venient variables can be taken on a suitably mathemati- 
cally well-founded formulation, there is no guarantee that 
numerical simulations will be stable, either around flat- 
space or the puncture data we are interested in evolving. 
Numerical tests will be discussed in detail elsewhere. 



VII. CONCLUSION 

Motivated partially by the success of symplectic inte- 
grators in many numerical applications we have studied 
Hamiltonian formulations of GR with live gauges. An es- 
sential property for any numerical application is that the 
underlying PDE problem is well-posed. Thus we have 
concentrated on well-posedness of the I(B)VP and con- 
sidered the level of hyperbolicity of the equations of mo- 
tion when the gauge choice is made at the Hamiltonian 
level. 

Several tools were developed to perform our analysis. 
We examined the relationship between the Hamiltonian 
energy and that required for symmetric hyperbolicity. In 
general we find that the Hamiltonian guarantees the ex- 
istence of a candidate symmetrizer, which is typically not 
positive definite in applications. However there are im- 
portant special cases in which Hamiltonian structure sig- 
nificantly simplifies hyperbolicity analysis. We have also 
translated the analysis in the literature on first order in 
time, second order in space systems to fully second order 
systems and now view the fully second order character- 
istic variables as the natural form whenever they can be 
constructed. 

In our first applications we demonstrated that every 
strongly hyperbolic Hamiltonian formulation of electro- 
magnetism is symmetric hyperbolic. We furthermore 
show that a generalization of electromagnetism, with 
Hamiltonian structure, the pure gauge system, is also 
always symmetric hyperbolic whenever it is strongly hy- 
perbolic. 

There are several interesting consequences of the 
Hamiltonian approach to gauge choice. Normally from 
the free-evolution PDEs point of view the gauge choice of 



GR (or electromagnetism) allows one to choose equations 
of motion for some quantities and the freedom to add 
combinations of the constraints to the equations of mo- 
tion. But in the Hamiltonian approach in GR for example 
the constraint addition is determined completely by the 
choice of the lapse and shift. This restricts the choice of 
formulation. There are choices of evolution equations for 
the lapse and shift that allow for a strongly or symmetric 
hyperbolic system without Hamiltonian structure that 
are forbidden when Hamiltonian structure is imposed, 
including the popular puncture gauge. We find in every 
example that the fully second order characteristic vari- 
ables of the Hamiltonian systems always have a special 
form. We do not know how generally this property holds, 
but there are certainly examples of strongly hyperbolic 
formulations without Hamiltonian structure that do not 
exhibit it. 

In our analysis of GR we find several families of 
strongly hyperbolic formulations and are able to choose 
the gauge parameters so that we get very close to the 
puncture gauge choice. Analysis of symmetric hyperbol- 
icity is rather more difficult, but for every strongly hy- 
perbolic formulation we are at least able to obtain par- 
tial results. We analyze several cases completely and find 
a Hamiltonian formulation that is symmetric hyperbolic 
with a small modification of the puncture gauge. 

The next practical step for applications is to investi- 
gate whether or not any of the formulations constructed 
here may be used successfully in numerical evolutions, 
either with or without symplectic integration. There are 
however delicate issues involved in the choice of evolved 
variables. 

An interesting question about GR is; what are the set 
of gauge choices that admit a well-posed I(B)VP? This 
question may be asked whether or not one restricts to 
formulations with Hamiltonian structure, and has still 
not been answered to our satisfaction in the general case. 
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Appendix A: Hyperbolicity of fully second order 
systems 

In this section we discuss the various definitions of hy- 
perbolicity for fully second order systems as summarized 
in Sect. QTI] 
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The key in the proofs is the following. If two matrices 



A ij = M- l [V ij - G'MF J ], 
B i = [F*M _1 + M~ 1 G l ]M 



are given then the matrix 

aV = 



o sf 

A pj B p 



is similar to the matrix A p i 3 : 



where 



= 

/ i — 

7? - 



ypi qp 
6 k 



51 
—MF j M 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 
(A6) 



This shows immediately that a fully second order sys- 
tem is strongly hyperbolic if and only if the fully second 
order principal symbol ([7]) has a complete set of eigen- 
vectors with real eigenvalues. 

Moreover one may easily prove that the notions of 
strong and symmetric hyperbolicity for a first order in 
time, second order in space system and an equivalent 
fully second order system agree. 

Now we show that symmetric hyperbolicity of the fully 
second order system is equivalent to the existence of a 
positive conserved quantity. 

Proof. Fully second order system is symmetric hyperbolic 
=> a positive conserved quantity exists: Without loss of 
generality we assume the following reduction to first or- 
der 



d t q = M w + F t diq + S q , 
d t w = V lj d l d j q + G l d t w + S w , 



(A7) 



which results in a first order in time, second order in 
space system with principal part matrix (|A4[) . We notice 
that in the principal part one can write 



dt 



diq 
d t q 



d 3 q 



T-VA',% ( d i q 



diq 
d t q 



(A8) 



It is then straightforward to check that if fP J is a sym- 
metrizer of the first order in time second order in space 
system then kH^Tv 1 defines a positive conserved 
quantity for the fully second order system. Positivity 



follows from Sylvester's law of inertia. To show conser- 
vation of the corresponding energy we notice that 

S{T^ i k H kl V^P i m s v S m = S^ l k H kl A p ^s p T 3 m S m 
= ■p s s S k H kl AP l i Sp S j T s s 
= -p s s (S k H kl A p l h p S j ^T s s 
= (Si-phH^A^hpTrSm^ 
= (S i Tt i kH kl TlW* j m 8 p S m ? i (A9) 

where Sj is an arbitrary spatial unit vector and 





1 



(A10) 



Using the techniques presented in [12| one can then show 
that (|A9p implies conservation of the energy 

e = \{drf ^)Tt^7p (g*Y (All) 



A positive conserved quantity exists =>■ Fully second 
order system is symmetric hyperbolic: This direction is 
obvious if one chooses the natural first order reduction 
with dtq — w. The symmctrizcr of the first order in time 
second order in space system is then the same as the 
symmetrizer of the fully second order system. □ 

Lemma 1. Symmetric hyperbolicity of the fully second 
order system J5J is equivalent to the existence of a second 
order symmetrizer and fluxes (Hi , (f> % , (f> lJ ) satisfying 



At* - 



(A12) 



= A^ k (<p l -H 1 B l )s l s jSk , (A13) 

for every spatial vector s,, and 

R x A lJ + B^HxBfl ~ B(«'t0j) + «p t <p - B^Hi 
ft - H t B j Ht 

(A14) 

hermitian positive definite. 

Proof. We proceed by demonstrating equivalence of 
(|A12[) - (|A14j) and the existence of a conserved positive 
quantity on the system. 

(j A12|) - (I A14|) =>■ conserved quantity exists: Consider 
the energy 



E = I d 3 x e, 



(A15) 



where is the matrix (|A14|) . The positivity of 
obviously implies positivity of E. It remains to show 
that E is conserved. 
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Computing a time derivative of the energy density e 
we get 

2d t e = d p (d t qU P d t q) (A16) 

+ d p (%U ipt #i% + dtq^H^diq) 
+ didpqU^H^ *d k q + d^H*A pk d k d p q, 

where H\ = tf-HiB*, Hf ] = ^t+gbt^i] -B^H^l 
The last two terms can be written as a divergence, too: 

didpq 1 A 1 ^ H k] 8 k q + d^H l 2 A pk d k d p q 

= didpqU^H^ ] d k q + d^H^A^dkdpq 

+ ^didpq^A^H^ -A k ^H^)d k q 

+ ^d iq \H l 2 A pk - H ( 2 k A^ l )d k d p q 



dp (d iq ^H^A pk) d k q 



di (d p q^(A lpJ< H kf - A k ^H^)d k q 



3 

+ \d k (d iq \W 2 A pk - H { 2 k A^)d p q) , (A17) 

where the last equality holds because of (|A13|) (for a 
tensor T Ijfe the equation T^') = is equivalent to the 
requirement T 1J ' fc SiSjS k — for every covector Sj). Hence 
_B is a conserved quantity. 

conserved quantity exists (|A12|) - (|A14|) : Assume 
that we have the conserved quantity (|A15|) with 

(A18 » 

and parametrize the fluxes of the energy by 

- - ( £S t ) ■ 

with (f)\ z = (p\ and 4>z^ k ^ = 03 lfe -'. Computing and com- 
paring dt£ and <9i$ l reveals that energy conservation im- 
plies 



2$ k)i = A jk ^H 2 ^, 
2<j>f> =A IJ] H 1 , 
2,q //;," • H 2 j B\ 

The solution of (IA20dl) is obvious: 
From (|A20bp and (fA20c) we get 



(A20a) 
(A20b) 
(A20c) 
(A20d) 

(A21) 



-H 2 l B 3 

4 1 



Since if y defines a symmetrizer we also have H l J ^ = 
fjg\ When we identify 

</>* = 24, (A23) 
<t> ij = 2cj) l 2 j] + WiB 3] - B^HxBfi 
this implies that (IA14[) is hermitian, and (IA12[) is obvi- 



Concerning (|A20a|) it was proven in [12, section IVB] 
that this equation together with 4> * — $ %k i has a solu- 
tion if and only if 



A {jk ^H 2 1 ^ = H 2 ^A 3k) . 



(A24) 



This equation is equivalent to (|A13I) . 

Finally, the positivity of the symmetrizer H lJ implies 
that the matrix (|A14|) is positive definite. □ 



Appendix B: Proof of the rank criterion 



In section ITVBl we described the rank criterion, a fea- 
ture that can be used as a necessary condition for sym- 
metric hyperbolicity. Here we discuss the proof of this 
criterion. 

We follow the notation of section ITVBl Thus, C l C M. 1 
is a set of formulation parameters, c, and we consider 
k x k principal part matrices A p i J (c) that depend con- 
tinuously on c. We search for symmetrizers of A p i J (c) in 
the set of hermitian k x fc-matrices, G, that is the image 
of a linear map G y : M™ — > R fcxfc , g i-> G'^(g). The re- 
quirement that G %3 (g) is a candidate symmetrizer for the 
formulation c defines linear equations on g of the form 

B(c)g = 0, (Bl) 

where B(c) E R mx " depends continuously on c. 
We prove the following 

Lemma 2. Let C l C C l be the set of formulations where 
the rank of B(c) is N and that do not possess symmetriz- 
ers in Q . Let further c be a formulation such that the 
intersection of every neighborhood U c C C l of c with C l is 
not empty, U c C\C l ^ 0. 

If there is a symmetrizer for c in Q then the rank of 
B{c) is smaller than N . 

For the proof we need the following 

Proposition 1. Let B,B G R mx " with \\B - B\\ < e 
and g G kerB. Let further be 



K(B) 



inf 



\Bv\ 



i>e(feerB)^\{0} \\V\\ 



(B2) 



A ij f i?i - 2cj) 2 ij] - 2<j)\B j + B^H^BK (A22) 



where V 1 - := {v G M" : (v, v) = 0, Vw G V}. 
Then we get 

inf Jg-h\\= mm_\\g-h\\<-^\\g\\ 

hekerB h&kerB K(B) 



(B3) 
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Proof. It is clear that there exists a. g £ kerB such that 
inf/refcers 11.9-^1 = \\9-9\\- It satisfies g - 5 € (kerB) x . 
If g = g nothing needs to be shown. Hence, we assume 
that g - g € {kerB) x \ {0}. The definition of K then 
implies 



K(B)\\g-g\\ < \\B(g - ~g)\\ = \\Bg\\ 
< ||B-B|||| fl ||<e|| 5 || 



\(B-B)g\\ 



(B4) 



□ 



Now we prove the above lemma. Since B : C l — > R mx ™ 
is continuous, there exists for every e > a S > such 
that \\B(c) - B(c)\\ < e for every c with \\c- c\\ < 5. We 
assume that g £ R ra satisfies B(c)g = 0. Then, according 
to the proposition there is for every c with \\c— c\\ < 5 a 
g £ ker(B(c)) that satisfies 



\\g-g < 



K(B(c)) 



9\\ 



(B5) 



Since U c r\C l ^ for every neighborhood U c of c we 
can choose c £ C l such that \\c ~ c\\ < S. Then we get 
for every g £ kerB(c) that G^{g) is not positive definite 
(because otherwise there is a symmetrizer for c £ C l ). 

Hence, for every e > and every candidate G^ (g) of 
c we find a candidate G lJ (g) of c that is not positive 
definite such that 



11.9 -.911 < 



K(B(c)) 



\9\\ 



(B6) 



On the other hand, since we find a symmetrizer for c 
in Q there exists a point g + £ kerB(c) such that G v (g+) 
is positive definite. 

We notice that G 1 - 7 is a continuous map (because it 
is linear) and that the map which assigns to a matrix 
its smallest eigenvalue is continuous, too. This implies 
that there is a neighborhood U+ C Q of g + such that the 
matrix G 1 ^ (h) is positive definite for every h £ U+ . 

We choose g+ £ ker(B(c)) such that 



< 



K(B(c)) 



1.9+1 



(B7) 



We know that g+ is not in U+. Therefore there exists a 
p > such that \\g + — g + \\ > p. 

Equation (|B7I) then implies that for every e > there 
is a formulation c £ C l with 



K(B{c)) < Ae, 



(B8) 



with the constant A = ||<jr_|_||/p. 

To show that this implies a jump in the rank of B we 
consider a sequence {c„}„ e N C C l such that 



lim c„ = 



(B9) 



We notice that lim^oo K{B(c n )) = 0. 

Then, for every c n there exists a g n £ (kerB(c n ))' L 
such that ||g„|| = 1 and K{B{c n )) = \\B(c n )g n \\ (the 



set {\\g\\ = 1} is compact). From the sequence {g n } we 
choose a convergent subsequence {h m } and get 

= lim K(B(c m ))h m = lim B(c m )h m 

m— foo m— >-oo 

= ( lim 5(c m )) ( lim h m ) = B(c)h, (BIO) 

\m— fcxD / \m— foo / 

with || ft || = 1. Analogously one can show that the limit 
of every convergent sequence {<?„ £ kerB(c n )} is in the 
kernel of B(c). 

It follows that the dimension of the kernel of B(c) is 
bigger than the dimension of the kernel of B(c n ), i.e. the 
rank of B(c) is smaller than N. □ 



Appendix C: Derivation of conditions for strong 
hyperbolicity 

In this appendix we discuss how conditions for strong 
hyperbolicity can be derived on the Hamiltonian formu- 
lations of GR introduced in section IVI Cl One finds that 
the special structure of the principal symbol which we as- 
sume in this article is very helpful for this purpose. We 
consider first order in time, second order in space formu- 
lations, but most steps can be directly carried forward to 
the fully second order case. 

The key will be to write the symbol in a simple stan- 
dard form where the conditions can be read off easily. 
This is meant to be a preliminary step for the construc- 
tion of positive definite symmetrizers in the next ap- 
pendix IdI 

We are interested in Hamiltonian formulations of GR 
with the Hamiltonian % given in (|80p . Thus, we are free 
to choose the seven formulation parameters Ci, . . . , G7. 
Yet, we replace C\, C2 and C3 by A 2 ^, Ay and \ 2 + respec- 
tively, using and (|88p . The given matrix expressions 
are valid in the Z4 variables (|111[) . But the analysis in 
the canonical variables is analogous. The only difference 
are factors of 2 at some places. 

As discussed in section fVI Cl for the hyperbolicity anal- 
ysis one can assume without loss of generality a = 1, 
f3 l = 0. I.e. the principal part matrix depends on the 
3-metric and the formulation parameters only. 

After changing the order of variables to 
(jij , a, Zi, Kij,Q, j3 l ) we find that the principal symbol, 
P s , has the block structure 



P s = 



X 
Y 



(Cl) 



To this matrix we apply similarity transformations to 
bring it to a form where the conditions for strong hy- 
perbolicity can be read off easily. First we consider the 
case where X is invertible. For our example that means 
A 2 . ^ 0. We get 



P s = T x l P s T x = 



1 

XY 



(C2) 
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where 



1 

o x- 1 



(C3) 



Since we assumed that the principal part matrix can 
be written in terms of the metric only one can write XY 
in block diagonal form. We define the orthogonal projec- 
tor, q, to the direction s and decompose the vectors and 
symmetric 2-tensors as 



V 1 = q\V A + s l V s 



q\Aq J B) - g 



V2 



(C4) 



where we use f sA := 8 i q j A T ii j\[2 = T sA /V2 and f ** := 
s/ZT^qtj — \plT qq instead of the usual decomposition 
with T sA and T", because it makes the transformation 
map orthogonal. We apply this decomposition to P s and 
get that it decomposes into a trace-free tensor, a vector 
and a scalar block, Pj, F , P v and PJ respectively. The 
corresponding submatrices of XY are 



(C5) 



XY\ / = q r ( x v ~ 2 ^ {x \^ 1 " 2Ce) ) ten) 



XY S 



where 
A 



A 
B C 



X% - (1 + C 5 )C 7 



(C7) 



C 7 /%/2 

(A+ - (1 + C 5 )(C 7 + A 2 _)) C 7 + \ 2 _ 



-, i A_ — C5C7 C7A^_/2 
( ~ 1 -2C 5 XI 



B 




C 7 +A 2 _+2A^-4A^ + 1 



(C8) 



2B n = 1 + C* 5 + 8C* 4 (1 + C 5 ) + 8C* 6 C 7 - 4A^ 

- 4(A| - 1)(C 7 + A^) + 3(C 7 + A^)A^, 
B12 - 8C 6 - 4(\ 2 V - 1) - (1 + C 5 )(C 7 + A 2 .) + A* . 

We do not get the vanishing upper right block in XYs 
when we start from the fully second order system. There- 
fore we prefer the first order in time system here. 

We see that the trace-free tensor block is always diag- 
onalizable, because q^qj) 1 — qijq kl /2 is the identity. The 
vector block has the structure 



(C9) 




\y 



It is hence diagonalizable with real eigenvalues if and only 
if 61 = 0, i.e. C 6 = (Ay - l)/2, and \ v > 0. 

Concerning the scalar part we use that the upper right 
block of XYs vanishes. A necessary condition for its 
diagonalizability is hence that A and C are diagonaliz- 
able. One finds that the eigenvalues of A and C are the 
same, namely ((89]) . Hence, A and C are diagonalizable 
if {\\ ~\ 2 _) 2 + C 2 C 2 + 2C b C 7 {\\ + X 2 _)^Q. 

When A and C have only one eigenvalue those matri- 
ces must be (as 2 x 2-matrices) already diagonal. This 
implies C5 = = C 7 . Moreover, since there is only one 
eigenvalue, the matrix B must vanish. One obtains (|96[) . 
The squared characteristic speeds are 

X v = 1/4(1 + 3A 2 .), v 2 +/ _=X 2 +/ _ = l + 2d, (CIO) 

and P s has real eigenvalues if A y > 0, A?_ > 0. However, 
the case X 2 _ — is not allowed here, because X would be 
singular. Hence, we get X 2 _ > 0. 

For v 2 ^ 7^ v 2 _ we can transform to a basis where A and 
C are diagonal. Let Ta and Tc diagonalize A and C 
respectively. We apply the transformation 



P s s = f^Plfs, 



(Gil) 



with 



Ts = <3iag(T A ,Tc,T A ,T c ) (C12) 
and get the scalar block PJ with the following structure 



/ 



XY, 




2 



_0 v 
B n B u v\ 
\B 21 B 22 v 2 _ 



(C13) 



Hence, to make Pg diagonalizable one needs that Bn = 
= B 22 - Those expressions are however very long so 
we do not present them. We discuss the solutions of the 
equations Bn = = B 22 in section fVl CI There we find 
three families of strongly hyperbolic formulations that we 
denote J1/2/3 ■ 

In the next section we prove symmetric hyperbolicity 
of the strongly hyperbolic formulations in the four pa- 
rameter family T\ (f9"3"|) which also satisfy C 7 = = C5. 
In that case we get 



Bn = {X 2 _-1- 8C 4 )/V4, 
B 22 = {3X\ - 8A y + 8C 6 + 5)/V2. 
Hence strong hyperbolicity implies 

X 2 V = (3X 2 + + l)/4, C 4 = (A^ - l)/8, 

v 2 + 



(C14) 



C 6 = 3(A 2 h - l)/8, 



XL > 0, 



a; >o. 



(C15) 



For X 2 ^ — the matrix XYs has vanishing eigenvalues. 
Hence, PJ is not diagonalizable. Thus, the case A 2 ^ = 
must be excluded. 
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Singular X : In the analysis above we needed the as- 
sumption that X is invertible. If X is not invertible 
(A?. = 0) then the diagonalizability of XY is only a nec- 
essary condition for the diagonalizability of P s (when P s 
is diagonalizablc then also (P 3 ) 2 , and hence XY and YX 
are). 

Indeed it turns out that the derived conditions are not 
sufficient in that case. However, having a complete set of 
eigenvectors for (P s ) 2 one obtains restrictions on the for- 
mulation parameters through the requirement that every 
eigenvector of (P s ) 2 must be an eigenvector of P s (de- 
tails can be found e.g. in fill ] ) . One finds a one parameter 
family of strongly hyperbolic formulations. But there the 
matrix M f|23[) is not invertible, i.e. we cannot transform 
to a fully second order free evolution system. Therefore 
those formulations are not of interest here. 



Appendix D: Construction of symmetrizers 

In this appendix a possible procedure for the construc- 
tion of a positive symmetrizer for formulations of GR is 
discussed. The difficulty is usually not the construction 
of general candidate symmetrizers but the proof of pos- 
itivity. We use a method that simplifies the latter task. 
The approach works as follows. 

First the principal symbol is transformed to a simple 
standard form using the calculations of the previous ap- 
pendix [Cj It is then possible to construct the set of all 
positive definite matrices, G, that symmetrize this trans- 
formed symbol. The next step is to make an ansatz for 
candidate symmetrizers and to compare the sub block of 
those ansatz candidates with the matrices in Q . If there 
is an overlap then the positivity of the matrices in Q may 
be used to simplify the positivity conditions of the ansatz 
candidate. 

It seems to be the most important task to find that 
kind of simplifications, because the reason why the proof 
of positivity is hard are the complicated expressions in 
the eigenvalues of the candidate symmetrizers. Starting 
from a symmetrizer of the symbol allows to combine some 
terms. 

We will show that in this way it is possible to prove 
symmetric hyperbolicity for the two parameter family of 
strongly hyperbolic formulations (|106l) (it is the special 
case of the four parameter family T\ (|Mf with C5 = = 
C 7 ). 

In section fVI CI we presented more strongly hyperbolic 
formulations, but we were not able to prove the existence 
of symmetrizers in the general case. The reason is that 
already the positivity conditions on the symmetrizers of 
the symbol, i.e. the matrices in Q, become complicated 
and we were not able to reduce them. 

Yet, for another two parameter family of strongly hy- 
perbolic formulations, (I107p . we were able to show that 
it is not symmetric hyperbolic, which means that in GR 
there are Hamiltonian formulations that are strongly but 
not symmetric hyperbolic. 



Decomposition of the ansatz candidate. In the con- 
struction of a symmetrizer we assume that we start from 
a reasonable ansatz candidate, G. The first step is a de- 
composition of its derivative indices into longitudinal and 
transverse part: 

G sB 
g ab 



G 11 = ( 



q-A 



G ss 

qAs 



(Dl) 



with G sst = G ss , G As ^ = G sA , G ABt = G BA . We 
can apply the same decomposition to the principal part 
matrix: 



SpA p ' 1 — 



A s 



A 



s B 



A S A S A s a b 



(D2) 



Due to the structure of A p i 3 we get s p q\A p ^ — 0, and 
it follows that G 1 - 7 is a candidate symmetrizer of A p j k if 
and only if for all spatial vectors s the matrix G SS A S S S is 
hcrmitian. 

Later we discuss the construction of a positive definite 
G ss . But given this matrix it is clear that G lJ is positive 
definite if and only if 



1 

£> q c 



G ss G sC 

qDs QCD 



G s 
G 



AB 



1 s^t 

q D B 


G As {G s 



IQSB 



with 



B A = _ G As( G ssy 



(D3) 



> 0, 



(D4) 



is. Hence, given a positive definite G ss we need to prove 
positivity of 



g ab _ G As (G ss )- L G 



IrtsB 



(D5) 



In our formulations of GR that means for the analysis 
of positivity the full candidate symmetrizer, a 40 x 40 
matrix, is reduced to two 20 x 20 matrices. 

Construction of a positive G ss . Now we want to con- 
struct a matrix G ss that symmetrizes the principal sym- 
bol P s = A s s s . One cannot expect that every sym- 
metrizer of P s is part of a symmetrizer for the full prin- 
cipal part matrix, i.e. the constructed G ss must be suffi- 
ciently general to make it possible to adjust parameters 
in the construction of the full symmetrizer later on. But 
it must not contain too many parameters, because this 
complicates the positivity analysis of the matrix (ID5[) . 

As a good compromise we found it reasonable to start 
from the ansatz that the full symmetrizer depends on the 
metric only, but is not restricted otherwise: 

G ijklmn = ( D g) 



/ s~iij kl mn 


kl 











kl m 
^16 


\ 


/-iji mn 
Lr 12 













^26 










riij , 
^33 km 


^34 k mn 


G35/C 













l - r 34m kl 


G44 kl mn 


G45 ki 













<-*35m 


G^45 mn 


G55 







\ G 16 


r<3 k 
°26 











r^k m 
U 66 


/ 
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with 



To simplify expressions we introduce new parameters: 



Gij kl ran 
11 



G 



ij kl 



U 22 



T 33fc rn 



G 



i kl m 
16 

Gi rn 
26 



34 A; mn 

Gi 
35k 

G44 

kl ran 
G45 kl 

Gk ra 
fifi 



Gl ~M ~M ~, mn 1 0/^2 „ij „,k(m nJ n)l 
n7 7 7 +zG n 7 J 7 v 7 ' 

+ 2Gl 1 (7 ^(fc 7 ^): '7 m,l + 7 i (' m 7 n W7 fei ) 

+ 2Gf 1 (7 ifc 7 j(m 7™ )i + 7 i 'y( m 7 ")' c 

+ A} 1 (7 fc [ 4 7 jl(m 7 n)z + y[*yK m y*) fn ) ) 
2Gl 27 l(fc 7 iW +G? 2 7 y 7 M , 
G 22 7 u ', 

2G^ 3 ^C) + Gi37' iJ 7 fem + A^M, 
2G ie7 v 7 y + ^ie7 7 , 

267 , 

G 34^( m 7n)fc + G^fcTmn, 

2G447fc( m 7 n ); + G| 4 7fc;7„ m , 

G 457 , 

Gl n ,km 
667 ■ 



Now, from the previous appendix [Q we know that for 
Ai 7^ there exists a matrix T and a diagonal matrix A 
such that 



1 

A 



(D7) 



We decompose the derivative indices of the ansatz can- 
didate according to (|D1[) and apply a congruence trans- 
formation using the matrix T: 



QSS 

qAs 



QSB 

G AB 
Tt 



(D8) 



q A C 



qCs qCD 



T 

q D B 



The condition that G 13 is a candidate symmetrizer of the 
full problem then becomes 



G ssps = (Qssps^ 



(D9) 



This is a linear equation on the parameters G c ab in G. It 
is guaranteed that the solution does not depend on the 
direction of s, because both the principal part matrix 
A p ^ and the ansatz candidate G u depend on the metric 
only. 

For the family (|106D the solutions of (|D9I) are five di- 
mensional. For the positivity analysis it is important to 
note that that G ss decomposes into 2x2 blocks. 



s v := 2G\ l + 



r 1 



2(3A 2 



1 



(D10) 



S44 



S22 



S12 



Gl 



^66 



33 2(3Ai + 1) ' 
4G? 1 -2Gi 3 + G| 4 -Gi 6 (l + 3A 2 + )- 1 



4G?i - 2G* 3 



Ai 

r 44 H 



G44 + 2G45 — ^66 (1 3A^_) 1 



2A 2 



where sy appears in the vector and S12, S22 as well as 
S44 in the scalar block. 

With those parameters and the observation that G ss 
decomposes into 2x2 blocks the positivity conditions on 
G ss can be reduced using computer algebra. One finds 
that for solutions of (|D9|) G ss is positive definite if and 
only if 



S44 < - 



S22 > - 



9A 2 S 2 2 



(A 2 _ S ? 2 



+ s 22 (G 2 X (4 + 9A^) - 2s v (l + 3X 2 + )) ) 



'12 



2(1 + 3A 2 _)(2G2 1 - Sy )' 



(Dll) 



and 



sv > 2G 2 n , 
S12 e K, 



S44 > 



2Gi\ 



XL > 0, 



G\i > 0, 
Ai > 0. 



To simplify that condition further we define new param- 
eters pi, P2, P3 as follows 



1 



Pi 



= tt (2s 22 (sy - 2G? X )(1 + 3A 2 ) - A 2 



a l : 



'12; ' 



P2 = N2/D2, p 3 = 
-N2 = S22(2G1 1 -.3 V )X 



3s 44 



2G? X 



- 1, 



(D12) 



D, = 



x (2G 2 X (4 + 9A+) + 9A 2 Ls44 - 4sy(l + 3A^)) 
+ {2G\ 1 ~ s v )2\ 2 _s\ 2 

4(G 2 !) 2 (4 + 9A 2 _)S 22 - §X\ S22S44SV 

+ 2G\ 1 {2\Ls\ 2 + 9X 2 + S22S44 - (4 + 9\ 2 + )s 2 2S V ) 



The positivity condition for G ss then becomes just 



Pi > 0, 
S12 e R, 



P2 > 0, 

x 2 _ > o, 



P3 > 0, G 2 ! > 0, 



\\ > 0. 



(D13) 



Positivity of the transverse block. Now we consider 
the remaining transverse block (|D5|) . It is easy to check 
that the transformation ()D8|) does not change this ma- 
trix. 

We use the solution of (|D9|) to write (|D5|) in terms of 
Pi, P2, P3, Gfi, S12, Ai, A^, A X1 and A\ 3 and show that 



22 



the parameters can be chosen to make it positive definite 



for every A 2 . > 0, A? > 0. 



First we notice that (ID5[) has the following structure 

iAB kl mn QAB kl 



rtAb 







gAH kl mn _ | G AB mn G AB q 







Gf 3 Bkm 



(D14) 



where the component tensors can be written in terms of 
q and s: 

GAB kl mn 1 ^AB kl mn , o 2 ~AB k{m n)l 

11 — 5ii<Z y y ~r z Snv y y 

+ 2gl 1 {q A ^ m q^ B q kl + q A(k q l)B q mn^ 
+ 2g\ l (q Ak q B{m q n » + q Al q B ^ m q^ k ) 
+ 2gl 1 {q Am q B{ - k q l ^ n + q An q B< - k q l ^ m ) 
+ g 6 11 q AB (s k s l q mn + q kl s m s n ) 
+ 2gl iq AB (s k s {m q n)l + s l s {m q n)k ) 
+ 2gl l {q A{ - m q^ B s k s l + q A( - k q l ^ B s m s n ) 
+ 2g 9 11 (q Ak q Bim s n h l + q Al q B(m S n h k ) 
■ 2g 1 1 1 (q Am q B(k S l h n + q An q B ( k s l h m ^ 



11 AB k I m n 

g n q s s s s , 



G 



AB kl 



(D15a) 

, - <7xW m + 2gf 2 q A ^ B + g 3 12 q AB S k s m , 

(D15b) 

(D15c) 



r AB - n 1 n AB 
U 22 ~ 9221 i 

r<ABkm _ 1 AB km , 2 „/lt Bra 

< '! 



T :s:S 



i 3 Am fcB , 4 m 
+ 533? 9 +533? s s ■ 



(D15d) 



When we assume that G22 is positive definite then we 
can transform this matrix without altering positivity to 



A;?* r*AB kl mn /~iAB ryABkm\ 
aia gOll , ^22 ) ^33 )' 



(D16) 



where 



GAB kl mn r^ABklmn r^ACkl(r^—l\ /~iDB mn 
11 •— Wl _ ^12 l U 22 )CD^\2 

(D17) 

The structure of Gf 1 Bklmn is of course ([DTKa) . too. Yet, 
we denote the scalar parameters g\ 1 instead of g\ x . 

Positivity of G AB . Concerning the matrix G AB we 
find 522 — ^-Pi/D> w ith the denominator 

D = 3G 2 u \l( Pl + s 2 12 ) P3 (l + pa) + 4pip 2 (3A^ + 1). 

(D18) 



We see that with the condition (|D13[) every term in N 
and D is positive. Hence, G AB is positive. 

Positivity of G AB km . If we write the matrix G AB k m 
in an orthonormal basis that contains s then we find the 
following eigenvalues 

{533, 933 ~ 533, 933 + 533, 033 + 2 ff33 + s| 3 }- ( D19 ) 



Written in terms of p%, P2, P3, G 2 l5 S12, A 2 ., X 2 + , A\ x and 
A^ 3 those eigenvalues have the form 

5 3 2 3 = 21(^33 -z 2 )(Al 3 -z 3 ), (D20a) 

533-533 = 43-^2, (D20b) 
933+533-^33 + ^3, (D20c) 
533 + 2 9 3 33 + 533 = «i04i ~ Z 5)(^33 ~ **)• ( D2 0d) 

The expressions for zx,...,z& are quite complicated func- 



tions of the parameters. But using (|D13[) and computer 
algebra one can show that that zi < 0, 24 < 0, z 2 < z 3 
and Z5 < Z3. Hence, we can choose A^ 3 such that G AB k m 
is positive definite, e.g. A\ 3 — (max(z2, Z5) + Zz)j2. 



Since A\ 3 does not appear in G AB or G AB kl '"" we im- 
pose conditions that do not restrict the parameter choices 
in the rest of the matrix. 

Positivity of G AB kl mn . The hardest part in the pos- 
itivity analysis is the matrix G ABklmn . Again the first 
step is to write it in an orthonormal basis that contains 
s. One finds that the resulting matrix decomposes into 
two 4x4 and two 2x2 blocks, where the 4x4 blocks 
are similar. The 2x2 blocks are 



n 7 n 10 

9 11 9n 

a 10 a 7 

9 11 9 11 



§1+ till + ?n 
5?i 



5ii 



til 
51 1 ? 



5?i 



(D21) 



Looking at the entries in the diagonal of the 4x4 blocks 
one finds complicated expressions only at two positions. 
Those entries are 

2(5n + til + 25h + 5ii + 5u) and g\\. 

Using a congruence transformation we decompose the 
4x4 blocks into 2x2 matrices such that the complicated 
expressions are in only one 2x2 matrix. The resulting 
matrices are 



9 I 9 2 ii + 5ii + 5n 5ii " 5u 
~' 5ii-5ii 5ii5ii+5ii 



1 / 2mn TO12 
tn + 25u V m i 2 TO22 



(D22a) 
(D22b) 



where 



mu = (tn) 2 + 5u(5n + 2^i) + 5ii(25n + ^9u + 5ii) 

-2((5 11 ) 2 + 25ii5ii+5ii(5ii+5ii)), 

mi2 = \/2(25ii5n " 2(5ii + 5ii)5ii + 5u(5ii + 5u)), 

m 22 = 5n5?i + 25i 1 i 1 5n " 2(5ii) 2 - (D23) 

To simplify the form of the eigenvalues of the 2x2 
blocks (ID21|) . (|D22al) we define another parameter, ^4, 
in terms of ^4^: 



Pi = 



G? f 2 + P2) (3A^ 3 (Pi + sf 2 ) + 4 P1 (1 + 3A 2 + )) 



A 1 



(D24) 
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We get the following eigenvalues 

4 (HGj^Xl + 3G 2 n (X 2 + - 1) P 4 - 4p|) 

3G 2 n - 2 P4 , 4 P4 , 6p 4 - 10p 4 - ^^}- (D25) 

Hence, together with (|D13p . the three matrices (ID21I) . 
(ID22a[) are positive definite if and only if the inequalities 

3 5 
p 4 > 0, PA < -G 2 U , Pi < -G 2 ^, 



P! > 0, p 2 > 0, 

si2 el, A^ > 



> 0, 



(D26) 



are satisfied. 

Now, to prove positivity of the remaining block (jD22b|) 
we use the fact that the only condition on s\2 is s\2 G M. 
We choose S12 such that (|D22b[) becomes diagonal, i.e. 
we solve the equation mi 2 =0 for s\2- We get sf 2 = 
—N12/D12 — pi with the numerator and denominator 

N 12 = 4(1 + 3X 2 + ) Pl p 2 x 

x [3(G 2 n ) 3 Xlp 3 + 3(G 2 n ) 2 X 2 + ((X 2 + - 1) P3 - 10) P4 
+ 4G? 1 (9-A»(p3-5))p2-24p|], 
r> 12 = 3A^p 3 [3(G? 1 ) 3 Al-P2P3 

+ 3(G? 1 ) 2 A^((A^ - l)p 2 p 3 - 10(1 +P2))pa 
+ 4G?i((9 + 5A* )(1 +P2) - X\p2Pz)pI 



24(l+p2)p2] 



(D27) 



The condition S12 € M can be written as s 2 2 > 0. It 
restricts the allowed range of P2 and p% : 



10p 3 > 



G n X + 



5G 2 -^ A 2 ^. 



\2 ,^2 \2 

6p 4 3G^ - 2p 4 p 4 



P2 



> 



6A^ 3 (5G^A^ - 6p4K3gji - 2 Pi ) Pi 
(4 + 3A^(4 + p 3 ))^2 



(D28) 



D 2 = 3(G 2 1 ) 3 X% P3 + 3{G 2 1 ) 2 X% ((A 2 + - l)p 3 - 10) p 4 
+ AG\ l (9-A 2 + ( P3 -5))p2_24^, 

but if these parameters are chosen sufficiently big then 
the inequalities are satisfied. 

If we replace S12 everywhere using (|D27[) then (ID22b[) 
has only one simple eigenvalue: 



2(5G 2 1 A 2 ;-6p4)(3G 2 1 -2p 4 )P4 
3(G 2 1 ) 2 A 2 ,+3G 2 1 (A 2 ;-l)p4-4p 2 



(D29) 



One can check that numerator and denominator are pos- 
itive if (|D26|) is satisfied. 

Hence, indeed for every A 2 . > 0, A+ > the parameters 
can be chosen such that G^\ B kl mn is positive definite. 

The calculations presented here were performed using 
the computer algebra system mathematica [27j with the 
package xTensor by Jose-Maria Martin-Garcfa [28j . 
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